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Recent Lattice Results on the Light Quark Masses 
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I discuss old and new determinations of the light quark masses using lattice QCD. Most lattice results 
using various approximations can be fit together in a simple picture which is consistent with lower values 
than conventionally supposed: in the quenched approximation for the strange quark mass in the MS scheme, 
ms(2 GeV) = 95(16) MeV, and for the average of the u and d quark masses, m((2 GeV) = 3.6(6) MeV. An 
estimate of the effects of light quark loops yields answers lower still: ms{2 GeV) in the range 54 - 92 MeV, and 
mi{2 GeV) in the range 2.1 - 3.5 MeV. 
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1. Introduction 

Among the most important applications of lat- 
tice QCD are the determinations of the funda- 
mental parameters of the standard model in the 
quark sector. Of these, one of the most important 
is the overall scale of the light quark masses. It is 
one of the most poorly known of the parameters 
of the standard model from prelattice methods. 
The Particle Data Group estimates a range of a 
factor of three in allowed values j^]: 

100 MeV < m,{l GeV) < 300 MeV, (1) 
5 MeV < md(l GeV) < 15 MeV, and (2) 
2 MeV < m„(l GeV) < 8 MeV. (3) 

(I will use rriq to represent the running quark mass 
in the MS scheme.) Global fits to standard model 
parameters are very sensitive to such large vari- 
ations. It is also one for which lattice methods 
provide the only systematically improvable deter- 
mination. This can be contrasted with the strong 
coupling constant a^, for example, for which it 
is possible to imagine going to higher and higher 
energy scattering experiments and extracting as 
with perturbation theory. Lattice quark mass ex- 
tractions are harder to do with solid error anal- 
ysis than as extractions, but they may be more 
important in the long nm. 

2. Prelattice Quark Mass Results 



2.1. Quark mass ratios 

The ratios of light quark masses can be investi- 
gated with some degree of reliability using chi- 
ral perturbation theory (xPT), which becomes 
asymptotically exact in the zero quark mass, zero 
energy limit of QCD . One combination of the 
light quark mass ratios is especially likely to be 
reliable, since it has been constructed to have van- 
ishing leading order corrections in xPT: 
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The other combination of quark mass ratios 
may be obtained at leading order in chiral per- 
turbation theory from the canonical prediction 
of xPT, that the meson mass squared is propor- 
tional to the quark mass: 
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The usefulness of this relation has been called 
into question by the discovery of a symmetry 
of the chiral Lagrangian which leaves physical 
predictions invariant under simultaneous trans- 
formations of the interactions, and of the quark 
masses 101 : 
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and cyclic in u,d,s. The effects on the pseu- 
doscalar spectrum of a nonzero u quark mass can 
be mimicked exactly by higher order interactions 
in the xPT Lagrangian and altered d and s quark 
masses. 

In light of this, Leutwyler has assembled a va- 
riety of arguments to test the size of the correc- 
tions to Eq. (H) and in particular the possibility 
of TO„ = Q . He concludes that the corrections 
are small and obtains 



ms/md = 18.9(8), 
mu/nid = 0.553(43). 
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2.2. Absolute value of the quark masses 

Chiral perturbation theory makes no state- 
ment about the third combination of light quark 
masses, the overall scale. A large range of results 
has been obtained with a variety of methods. See, 
for example, a compilation of results for from 
Ref. shown in Table [l[ The most reliable of 
these prelattice results are perhaps those using 
QCD sum rules, but even here, the systematic 
errors of the method are hard to pin down. It is 
here that there is the largest spread in existing 
results and that lattice methods probably have 
the most important role to play. 



Table 1 

Compilation of results for TOs from Ref. |l] . If de- 
fined, the renormalization scheme is rns{l GeV). 



nis (MeV) 


method 


194 (4) 


quark model 


118 


sum rules 


175 (55) 


sum rules 


> 300 


sum rules 


112 (66) 


xPT -I- estimate of (qq) 


378 (220) 


XPT -I- estimate of (qqqq) 


150 


strange baryon splittings 


135 


SU(6) 



3. Lattice Quark Mass Determinations 

Lattice determinations of standard model pa- 
rameters require: 1) fixing the bare lattice pa- 
rameters from physics, and 2) obtaining the MS 



parameters from these with short-distance match- 
ing calculations. Spin averaged splittings in the ip 
and T systems are convenient quantities to set the 
lattice spacing ||^ . The mesons are small and easy 
to understand because the quarks are nonrela- 
tivistic. Light pseudoscalar meson masses are the 
most convenient quantities to fix the light quark 
masses. Chiral symmetry makes them very sen- 
sitive to the quark masses, the mesons are small, 
the correlators have good statistics and are easy 
to fit over long time separations. The correlators 
have simple behavior as the quark mass is var- 
ied toward zero, in contrast with unstable vector 
mesons. Other quantities must give the same re- 
sults as the approximations are removed. 

The second piece is the determination of the 
parameters of the MS Lagrangian by matching 
perturbative, dimensionally regularized short dis- 
tance amplitudes to their lattice counterparts. It 
is desirable to do the lattice part of such calcula- 
tions nonperturbatively as much as possible, to 
test for the presence of nonperturbative short- 
distance effects and possible poor convergence 
of perturbation theory. Such nonperturbative 
short-distance calculations are harder to design 
for quark masses than for the strong coupling con- 
stant. Such nonperturbative short distance anal- 
ysis of quark mass extractions is currently less ad- 
vanced than the analogous investigations for the 
strong coupling constant as- 

The perturbative expression giving the MS 
mass from the lattice bare mass toq may be writ- 
ten 
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The mean-field-improved bare mass to is given 
by TO = mo/(l — ]^3^) in perturbation theory, 
and rh = toq / \/ {Up)mc nonperturbatively, if 
the expectation value of the plaquette, (J7p), is 
used to define mean field improvement ||6[ . 70 = 
l/(27r^) is the leading quark mass anomalous di- 
mension. This coefficients in this expression are 
well behaved for the Wilson action and the 0{a) 
improved action of Sheikholeslami and Wohlert 
(SW) § (for which Cm = 1-67 g and 4.72 § re- 
spectively). For staggered fermions. Cm is 132.9 
leading to renormalization factors of 50-100% 
which are not explained by large tadpole graphs. 
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Figure 1. Old lattice quark mass results were 
reviewed in Ref. Quenched staggered re- 

sults (dark squares) show good cutoff indepen- 
dence, but suffer from huge perturbative correc- 
tions. Quenched Wilson results (circles) had well- 
converged perturbation theory, but large cutoff 
dependence. The white squares show some un- 
quenched staggered results. 



The status of quark mass extractions as of 
a few years ago was reviewed in Refs. JTo| , pl| . 
Some data tabulated in Ref. [|oj are shown in 
Fig. ^ Unimproved perturbation theory was 
used, and the lattice spacing was determined with 
the p mass. The quenched results for staggered 
fermions (dark squares) are relatively indepen- 
dent of the lattice spacing, but Eq. |l^ is untrust- 
worthy because of the very large quantum cor- 
rection. The quenched Wilson results (diamonds) 
have much better controlled perturbation theory, 
but are much more lattice spacing dependent. If 
the results are extrapolated in a, they approach 
the staggered results more closely. However, the 
magnitude and origin of remaining lattice spacing 
errors is then unknown. Improved actions must 
be used to investigate this. 



4. Recent Quenched Lattice Results 

4.1. Unimproved Wilson fermion results. 

New unimproved results have been presented 
in the JLQCD contribution to this volume by 
T. Yoshie iQ. They include a data point at 
(3 — 6.3, corresponding to a smaller lattice spac- 
ing (a = (3.29 GeV)^^) than previous results. 
Fig. H shows their Wilson fermion results su- 
perimposed on a subset of the world data. Mp 
has been used to set the lattice spacing. Im- 
proved perturbative theory has been used in the 
renormalizations, which makes the Wilson results 
slightly higher than in the older analyses, and the 
staggered results much higher, about 50%. The 
new JLQCD Wilson results are statistically con- 
sistent with previous results, and so can be com- 
bined with them for a consistent extrapolation of 
the leading error. Taken by themselves, however, 
they appear to have a somewhat smaller a depen- 
dence then other results. 

Wilson data were also presented by Fermilab 
(upper points in Fig. ||) |l^,|lj]. These points 
also lie slightly below and have a smaller slope 
than most of the world data in Fig. ||. This is 
due to the fact that the Fermilab data use the 
spin averaged IP-IS splitting in the system to 
determine the lattice spacing, which is expected 
to have a small 0{a) error. Determining a with 
the p mass is equivalent to calculating mi /nip. 
Since the ipapi^Fp^^ip C'(a) correction operator is 
expected to push spin partners like the tt and the 
p apart, 0{a) errors in both masses contribute 
to the a dependence. When the IP-IS splitting 
is substituted for the p, the slope ought to be 
reduced, as observed. If p mass lattice spacings 
are substituted into the Fermilab data, the results 
line up with the upper points in Fig. ^. 

Gupta and Bhattacharya have performed a lin- 
ear extrapolation on the world Wilson data and 
obtain m,(2 GeV) = 3.1(3) MeV ||. They use 
slightly different analysis methods from those of 
Ref. and obtain a result below that which 
would be obtained from Fig. |^, which is about 
3.8 MeV. However, the authors say that their er- 
ror estimate is a "guess chosen to reflect the var- 
ious uncertainties discussed . . . and is not based 
on systematic analysis" . 
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Figure 2. Recent unimproved Wilson fermion re- 
sults for TO/ (2 GeV) from JLQCD compared with 
previous world data (upper points) p^ . Lattice 
spacing from Mp. JLQCD staggered fermion re- 
sults are the lower *'s. 



JLQCD also presented new results for stag- 
gered fermions which left the existing situation 
unchanged. The results are almost lattice inde- 
pendent for a < (1.5 GeV)^^, but it is difficult to 
estimate perturbative uncertainties since the per- 
turbative renormalization factor is so large. (At 
/3 = 6.4, they use Z = 1.69, a 69% correction of 
which most is unexplained by mean field improve- 
ment.) 

4.2. Improved quenched Wilson fermion 
results 

To investigate the size of the remaining cut- 
off errors in the Wilson action {0{a^), 0{aas), 
0{a'^), etc.), it is necessary to remove the leading 
error with an improved action. In the SW action, 
the leading 0{a) error is removed with the addi- 
tion of the operator ipa^i^Fpi^ip. There are several 
gauge links in the operator, so tadpole improve- 
ment predicts a rather large correction, of order 
50%. 

New results with the tadpole improved SW ac- 
tion were presented by Fermilab [|l3| . The lattice 
spacing dependence observed in the Wilson points 
is substantially reduced, but not completely elim- 
inated. If this arises from a residual 0{a) er- 
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Figure 3. Fermilab light quark masses for Wil- 
son action (circles) and the tadpole improved SW 
action (triangles) jl3|. Lattice spacing from the 
IP-IS splitting in the V' system. 



ror, the results extrapolate down to rather near 
the staggered results. (A recent nonperturbative 
determination of the coefficient of the SW im- 
provement operator has indicated the possibility 
of even larger corrections than the sizable ones 
indicated by perturbation theory and mean field 
improvement fl^].) If it arises from higher or- 
der effects in a, the answer is close to the exist- 
ing small lattice spacing points, and the discrep- 
ancy with staggered results must be attributed to 
poorly behaved staggered perturbation theory. 

Since a downward trend in a is still present in 
the tadpole improved data, we take the improved 
result at the finest lattice spacing as an upper 
bound in the quenched approximation and a lin- 
ear extrapolation as a lower bound. Adding per- 
turbative uncertainties linearly and other uncer- 
tainties in quadrature gives the quenched result: 



to; (2 GeV) = 
TO, (2 GeV) = 



3.6(6) MeV, 
95(16) MeV. 



(11) 
(12) 



Another determination of the strange quark 
mass with an 0{a) improved action has been re- 
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Figure 4. (aM^)^ vs. aruq for staggered fermions 
at /3 = 6.4 on a 40^ x 96 lattice @. Deviations 
from linearity are a few % or less. 



ported, in Ref. This determination used a 

tree-level, rather than a mean-field improved co- 
efficient for the improvement operator. They ob- 
tained m, (2 GeV) = 128 ± 18 MeV. They did not 
attempt to correct for the effects of the remaining 
lattice spacing dependence or the effects of the 
quenched approximation. Much of the discrep- 
ancy with the Fermilab results arises from fact 
that we have used much larger SW improvement 
coefficients, and make an allowance for the fact 
that we continue to find significant cut-off depen- 
dence even so. 

5. Corrections to Quark Mass Ratios 

5.1. Nonlinearities in niq vs. A/^ 

There is very little lattice evidence for large 
deviations from linearity in the quark masses in 
Eq. (H) in either quenched or unquenched calcu- 
lations. In the most accurate data, the question 
often seems to be why the predictions of chiral 
symmetry appear linear at the few % level up 
to such large values of the quark masses. (See 
Fig. IJ. The largest pion mass is around 850 MeV, 
and nonlinearity in (aM^)^ vs. arUq is only a 
few %.) However, searching for such deviations 



is tricky, especially in the quenched approxima- 
tion where quenched chiral logarithms may add 
spurious nonlinearities at small rrig fl9|] . Exist- 
ing unquenched calculations have not yet exam- 
ined carefully the case of broken flavor SU(3), 
nis > mi, so not all of the higher order operators 
of the chiral Lagrangian have yet been tested. 

5.2. Electromagnetic effects 

It is in principle simple to determine the con- 
tributions of electromagnetism to hadron masses 
using numerical simulations with link matrices 
which are products of SU(3) and U(l) matrices. 
The U(l) phases of electromagnetism are cheap 
to generate because to leading order in aem (all 
that is required for practical purposes), they can 
be obtained from Fourier transforms of Gaussian 
fluctuations in momentum space. To see the ef- 
fects of electromagnetism clearly above fluctua- 
tions in the SU(3) fleld, one wants to use values 
of aem which are larger than the physical value. 
It then remains only to show that electromagnetic 
effects are still linear in this region. That this in- 
deed holds has been shown recently by Duncan, 
Eichten, and Thackcr ||2^. They have performed 
a prototype calculation of the 7r+ — 7r° splitting at 
P = 5.7 and with it have obtained mu/rrid = 0.51, 
in agreement (so far) with Eq. ^. 

5.3. Can m„ = 0? 

The most interesting application of these cal- 
culations is the settling of the question of whether 
m„ = in the real world. This possibility is fer- 
vently desired in spite of all evidence to the con- 
trary because of its neat solution to the strong 
CP problem. Rcf. § showed that allowing cor- 
rections as large as 30% to xPT equations such 
as 1^ made m„ — compatible with meson mass 
data. (See Fig. ^, however.) The possibility of 
large instanton-induced flavor mixing effects has 
been proposed as a mechanism to generate such 
corrections in QCD ||2^ . Not all the required lat- 
tice calculations have been done, but so far, the 
lattice evidence is against m„ = 0. 

6. The Quenched Approximation 

As the lattice spacing is reduced, while keeping 
hadronic physics flxed, the lattice bare couplings 



evolve according to their anomalous dimensions. 
In the quenched approximation, these anomalous 
dimensions are slightly wrong, due to the absence 
of light quark loops. The strong coupling con- 
stant evolves according to the zero-quark (3 func- 
tion coefficient /3q*''' = 11 rather than the correct 
three-quark coefficient p'^^ = 11 — 2/'i^nf = 9. 
Asymptotically, ctsij^ jo) in the quenched approx- 
imation is expected to be too small by a factor of 
about 9/11. 

Since the short-distance quark mass evolu- 
tion is given by d In m(q)/(i In (g) = — 7oas/(47r), 
where 70 = 8, this implies that the quark mass 
evolves too slowly in the quenched approxima- 
tion, and at small lattice spacings is larger than 
in real life. At high energies, running mass evo- 
lution is given by 



To leading log accuracy, therefore, the effect of 
the absence of quark loops due to perturbative 
effects on the evolution of the running mass from 
the strong coupling region (where « 1) to the 
high energy region can be approximated by p2[ 



m(7r/a)| 



qu. 



m(7r/a)[unqu. 



a{n/ay^(^/^o~^/€') (14) 
1.1 to 1.2, (15) 



for a{n/a) w 1/4 to 1/8. 

There are also quenching effects arising from 
the nonperturbative region. Unlike the case of 
quarkonium systems, for pions there is no argu- 
ment that these should be smaller than the per- 
turbative effects. The above expression can't be 
taken as a correction factor, only as an indication 
of the direction and order of magnitude of effects 
to be expected. Nonperturbative calculations are 
required to investigate quenching effects quanti- 
tatively. 

Unquenched results for Wilson fermions appear 
complicated and hard to interpret. They differ by 
as much as a factor of two from quenched results, 
and do not seem to reproduce the lattice spac- 
ing dependence of the quenched results as would 
have been expected (see Ref. [0). Since most 
other unquenched Wilson results, such as those 



in thermodynamics, are also hard to understand 
compared to staggered results, I will not attempt 
to fit unquenched Wilson fermion results into my 
general picture. 

Some unquenched staggered results summa- 
rized in Ref. [0 are shown in Fig. |^ (white 
squares) along with the quenched results. The 
unquenched results indeed lie below the quenched 
results by roughly the expected amount, and I 
will use them to estimate the effects of quench- 
ing. The large corrections in staggered fermion 
mass renormalization cancel out in the ratio of the 
quenched and unquenched determinations, mak- 
ing this an useful quantity to examine. To min- 
imize effects due to differences in analysis meth- 
ods, we estimate the ratio from the results of a 
single group, at similar volumes and lattice spac- 
ings (about 0.4 GeV"^) pSlpi], and obtain 



(13) m;(1.0 GeV)n,^o 



2.61(9) 



m/(1.0 GeV)nf=2 2.16(10) 
= 1.21(7) 



(16) 
(17) 



Since there are, in fact, three flavors of light 
quarks in the world and not two, I will this ratio 
as a lower bound on the actual ratio and use the 
square (corresponding to four light quarks) as an 
upper bound. 

7. Synthesis 

Lattice determinations of light quark masses 
are more difficult than the analogous determina- 
tions of as ■ Pion masses have worse a dependence 
than the quarkonium splittings. Finding nonper- 
turbative ways of eliminating the perturbative re- 
lations with the bare lattice quark mass requires 
more work. 

Nevertheless, with a few exceptions, most lat- 
tice quark mass extractions are consistent with a 
reasonably simple picture. I ignore results with 
very large a or very small volume. Unquenched 
Wilson results, which are also hard to interpret in 
most other quantities, do not seem to make sense 
compared with quenched Wilson results. Of the 
remaining results, the larger magnitude and lat- 
tice spacing dependence of Wilson results com- 
pared with staggered results is greatly reduced 
with 0{a) improved actions. The magnitudes are 



reduced more with tadpole improved SW correc- 
tion terms than with tree-level correction terms. 
If the remaining discrepancy arises mainly from 
residual 0{a) effects in the improved action, the 
true quenched answer lies close to the staggered 
result, nil a little over 3 MeV. If it arises mainly 
from higher order corrections in the staggered 
fermion lattice-M5' mass conversion (where the 
leading correction is 50-100%), the true answer 
lies closer to the improved result, m; ~ 4 MeV. 

Unquenched staggered results lie somewhat be- 
low quenched staggered results, but by an amount 
which is reasonable. Taking the ratio from Sec. |^ 
and using it to make a correction on our quenched 
result we obtain 

• ms{2 GeV) in the range 54 - 92 MeV, 

• m;(2 GeV) in the range 2.1 - 3.5 MeV, 

for the AIS masses renormalized at 2 GeV. The 
uncertainties most in need of further study are 
those associated with lattice spacing dependence 
and the quenched approximation. 
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